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Abstract 

One of the most novel applications of the Stirling cycle is in the free piston configuration that was 
initially designed by W. Beale. In free piston Stirling engines (FPSEs), there are no mechanical linkages 
coupling the pistons or displacers, the motions of the reciprocating components follow the working gas 
pressure variations. Fillipo de Monte and G. Benvenuto have recently proposed a linearization technique of 
the dynamic balance equations. The aim of this paper is to predict the thermodynamic conditions for stable 
operation of FPSEs and their modeling. 

The equations of the angular velocity are solved analytically in terms of the working gas mass and the 
displacer-piston phase angle of the machine. Using the criterion of stable engine cyclic steady operation, a 
mathematically rigorous form is obtained for the main parameters of the engine. 

Furthermore, for simplicity reasons, thermodynamic magnitudes are obtained using the Schmidt analysis 
(isothermal model). 

© 2003 Elsevier Ltd. All rights reserved. 
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1. Introduction 

One of the most novel applications of the Stirling cycle is in the free piston configuration that 
initially was designed by W. Beale at Ohio University in the late 1960s [1,2]. This configuration is 
the one that holds the most immediate promise. Free piston Stirling engines (FPSEs) have no 
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Nomenclature 

A cross-sectional area for moving element (cm 2 ) 

C c compression clearance space (mm) 

C Hdc displacer gas spring damping (N s/m) 

C Hpc piston gas spring damping (N s/m) 

C pc load damping (N s/m) 

D damping coefficient (s ' ) 

E e expansion clearance space (mm) 

F force (N) 

K c casing stiffness (kg/s 2 ) 

M mass (g) 

rij eigenvectors 

P power (kW) 

p pressure (bar) 

Q heat transfer (J/cycle) 

R working gas (helium) constant (J/kgK) 

r displacer-piston amplitude ratio, X A /X V 

r, eigenvalues 

S stiffness coefficient (s~ 2 ) 

T temperature (K) 

t time (s) 

V volume (cm 3 ) 

v velocity (m/s) 

W work done (J/cycle) 

x displacement (m) 

X displacement amplitude (m) 

Greek symbols 
A difference 

<j) piston-displacer phase angle (degree) 

to angular velocity (rad/s) 

Subscripts 

B bounce space 

C compression space 

CL clearance space 

c casing 

d displacer 

E expansion space 

g gas spring 

H heater 

K cooler 
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mean average 
p piston 

R regenerator 

Superscripts 

first order derivative with respect to time 
second order derivative with respect to time 


kinematic mechanism coupling the reciprocating elements to each other or to a common rotating 
shaft. Instead, the elements move entirely in response to the gas or other spring forces acting upon 
them. There is a surprising diversity in the various types of FPSEs [3]. The most common 
arrangement used for useful applications is a single cylinder (b-type). 

FPSEs can be used as thermal engines as well as cryocoolers. As prime movers, they use dif¬ 
ferent types of fuel: gas, liquid and solid in some cases. Nuclear [4,5] or solar energy [6] may 
be utilized also as a power source. This machine usually is connected with a linear alternator (LA) 
to produce electric power [7], 

There are two distinct preferred modes of operation (single cylinder, b-type): (i) Large piston 
motions relative to casing motions, when the power is mainly from the piston. In this case, the 
mass of the piston is lighter than that of the casing, and the piston amplitude is greater than the 
casing amplitude. In this group belongs the RE-1000, which was designed by Sunpower Inc. some 
years ago. (ii) The other mode is large casing motions relative to piston motions, when the power 
is mainly derived from the casing, and its mass is lighter than that of the piston. Respectively, the 
casing amplitude is greater than the piston amplitude. Such type of machine is the M-100, which 
was also designed by Sunpower Inc. 

The main advantages of FPSEs relative to kinematic Stirling engines are: (i) simpler mechanical 
design; (ii) the possibility to avoid any side force on the piston and also, by using flexure bearings, 
the possibility to eliminate all rubbing parts, reducing engine wear [8]; (iii) a high energy con¬ 
version (thermal to mechanical) efficiency; (iv) easy starting; (v) high performance; (vi) very long 
life; and (vii) low cost. The principal disadvantage of the FPSE is simply the lack of a rotating 
shaft. So many machines in engineering involve rotating shafts that their absence in FPSEs is 
sometimes regarded as a disadvantage, although in practice, the disadvantage is not so great as it 
appears. Many systems involving pumps, compressors and other machinery include kinematic 
mechanisms whereby the rotating motion of driven shafts is converted to reciprocating motion. 
All these can be directly driven by FPSEs and so can a reciprocating pump providing high 
pressure hydraulic fluid to drive hydraulic motors (e.g. in vehicles) [3], 

The correct choice of piston, displacer and casing masses, gas springs and other components are 
necessary for stable dynamic behavior. On all FPSEs, the displacer is typically much lighter than 
the piston. For positive power, the piston and casing motion are required to lag behind that of the 
displacer. FPSEs have load dependent frequencies. However, it will be shown that in some de¬ 
signs, with prudent selection of geometry, it is possible to reduce the variation of frequency with 
load to a negligible level [9], 

In most FPSEs, the gas flow is predominantly turbulent in the cooler and heater and laminar in 
the regenerator. The pressure drop in the machine is the summation of the pressure drops in the 
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cooler, regenerator and heater. However, the variation of density with pressure is neglected, since 
in most high pressure FPSEs, the pressure swing is usually a small fraction of the charge pressure, 
typically less than 15%. It is, therefore, assumed that density is only a function of local temper¬ 
ature, which, in the isothermal case, is constant. 

Optimizations of the dynamic behavior of FPSEs are not a simple matter because of the lack of 
mechanical linkages that are able to fix strokes (displacer and piston) and phase angle for the 
moving elements (displacer and piston). 

Because the laws of motion of the moving elements (displacer and piston) of FPSEs are not 
known presumptively, as in the case of kinematic Stirling engines, they have to be calculated by 
means of an accurate dynamic analysis. If the equations of motion are linearized (linear dynamic 
analysis), there will be a system of linear and homogeneous differential equations with constant 
coefficients. These equations may be solved analytically in the time domain. Solution of the 
polynomial characteristic equation yields the condition for periodic steady operation of the 
FPSEs [10]. 

An attempt has been by Rogdakis et al. [11], in which the equations of motion have been solved 
analytically in terms of stiffness and damping coefficients. 


2. Analysis 

2.1. Generalized analysis 


In order to describe the action of FPSEs, it is useful to assume harmonic motion of the 
components, which is often nearly true. The mass-spring-damper system in Fig. 1 illustrates 
the statements made above. The damping force represents work done by the mass, and as result 
of this, the applied force F must have a component normal to the displacing surface (parallel 
to the velocity) pointing up. These three forces (the applied force, the force due to the spring 
effect and the force due to the damping effect) are collinear, as it seems in the following equation 
[ 12 ]. 

Newton’s second law yields: 

F T F pnng + F damping — M • X ( 1 ) 


In a simple FPSE, there are three oscillating components, the piston, displacer and casing. If S tJ 
and D,j account for the influence on the z-component motions by virtue of a spring and damping, 
respectively, coupling to the motions of the /-component, whereas S n or D„ are a spring or 
damping effect unique to the /-component, respectively, the forms of the dynamic equations of the 
component parts (per unit element mass) are generalized as follows: 


Vp — .S'ppVp -(- *Sp d v d T .SpcXc -f /JppXp A E)p d x d T D pc Xq 


( 2 ) 


x d AipXp T AidTi T Alc-tc T D dp x p T Fd d x d T F dc x c 


( 3 ) 


A-pXp ~\- *S cd Vd “t“ 5' cc Xc \ /-TpXp “t“ ZTdX(| ~\~ D CC X C 


(4) 



te a Dove equation system may atso oe written 

A'pp Ap d »Apc Vp A^pp 

■Alp 5dd -A| c ' Vp ~t~ /-Alp /A 

-Ap -A\i -Ax- d\ -D cp / A 

the case of a FPSE of the b-type, the equation 
ritten: 

ir the piston: 

/V/pAp (pc Pb )'4 p ( A | [p c C pc )(Ap +x c ) 

he displacer: 

M d x d = A d Ap + A r ip c - Pd) ~ C Hdc (x d + x c ) 
he casing: 

MJA = Aa(dh 4- Ad — On) 4- Ao(n„ — da) — (C „„ 
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where the expansion space pressure is taken to be given by 

PE = Ap + pc (9) 

Because of the predominantly piston motion mode, the casing motion is ignored, and the 
describing equations are: 

For the piston: 

MpXp = (Pc Pb)^p (^Hpc T Cp C )Xp (10) 

For the displacer: 

M d x d = A d Ap + A R (p c - - C Hdc x d (11) 

The operation of the engine is described by the isothermal analysis (Schmidt analysis) [13], 
Using linear dynamics theory [9] that adequately describes the behavior of the FPSE, the spring 
and damping coefficients of the above equations as functions of the working gas mass M, angular 
velocity to and piston-displacer phase angle cf>, which is the phase angle between the motion of the 
piston and the motion of the displacer, are: 

5pp(M) = “ M ^ mean(M) {t^s + v b ) (12) 


S P d(M) 


Sdp(M) 


M * s PmUM) 



A pA R/^ m ean (M ) 

m d t k s 



S dd (M) 


A R/?mean (M ) 
Mr, 


Ap (Ap — A r ) A r 
7k S T k S 


°pp — 



M p 


D pd = 0 


D dp (M, to, </>) 


C P (M, co, </>) 
My, 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


C d (M,co, </>) — Cn d( 

Mr, 


ApCc i V R , F R ln[7k/r K ] i fk , A d E E 


D dd (M,co, <f>) = 

where 
5 = 


T + f«r + - 


+ ^ + - T 
r h r H 


7k -7k 

Now, the equation system (5) can be represented as follows: 
x = A • x 


(19) 


( 20 ) 


(21) 
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where 


x = 


and 


x = 


A = 


x d 

JCp 

x d 

' D pp 
-Ddp 

1 

0 


Xd 


x d 


-Dpd 

-Odd 

0 

1 


( 22 ) 


^pp 

5 dp 

0 

0 


5 pd 

Sdd 

0 

0 


(23) 


(24) 


The solution of the above system of differential equations is based on the roots of the char¬ 
acteristic polynomial: 

G(r) = r 4 +a- r 3 + p-r 2 +y-r + S (25) 

where 

a = D dd (M , to, </>) + D pp (26) 


P = D Pp • D dd (M, co, </>) - D pd ■ D dp (M , co, </>) - S dd {M) - 5 PP (M) 


(27) 


7 = D dp {M, co, </,) ■ S pd (M) + S dp {M) ■ D pd - D dd {M , co, </>) • S PP {M) - S dd (M ) • (28) 

and 

<5 = S PP (M) ■ S dd (M) - S pd (M ) • S dp (M) (29) 

The eigenvalues of A are r x , r 2 , r 3 and r 4 , with corresponding eigenvectors n i? n 2 , n 3 and n 4 . 
Therefore, the general solution of Eq. (20) is 

x(t) = ci • e r, t • ni + C 2 ■ e r2 > • n 2 + c 3 • e ryl • n 3 + c 4 • e 1 ' 4 ' 1 ■ n 4 (30) 

Because the vectors n,, n 2 , n 3 and n 4 are linearly independent, it is possible to find a unique set 
of constants c\, c 2 , c 3 and c 4 to satisfy any initial condition. The initial values x(0) imply that 

x(0) = ci • nj + c 2 • n 2 + c 3 • n 3 + c 4 • n 4 (31) 

According to the form of these roots, we have the following cases: 


(a) Roots with positive real parts: The exponential functions tend to infinity with time, and there 
is not any oscillation. 

(b) Roots with only negative real part : The equations of motion of the piston and displacer tend 
to zero. If there are two complex conjugate roots, there are convergent oscillations. 
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(c) Two roots with negative real parts and two roots with positive real park. x p {t) and x d {t) increase 
with the time. If the roots with positive real parts are complex conjugates, there is an oscilla¬ 
tion for each component, although they are divergent. 

(d) Two imaginary roots and two roots with a negative real part : Steady oscillations. 


The main conclusion is that the machine works at a cyclic steady state when the characteristic 
polynomial has two imaginary roots and two roots with a negative real part. The absolute values 
of this real part must be very high to ensure a quick drop of the transient and, therefore, a quick 
arrival of the stationary oscillations for the reciprocating parts. 

For secure stable operation [10,14] of the engine, the roots of the characteristic polynomial 
(Eq. (25)) should be 

r\ = k ■ i, r 2 = —k ■ i, r 3 = / and r 4 = m (32) 

According to the de Moivre form, these should effect 

(r — k ■ i)(r + k ■ i)(r — l)(r — m) = 0 =4> r 4 — {m + l)r 3 + {ml + k 2 )r 2 — k 2 {m + l)r + k 2 ml = 0 

(33) 

The above polynomial (Eq. (33)) must be equal to the characteristic polynomial G{r), (Eq. 
(25)). From this requirement, the following equations are obtained: 


—{m + /) = a (34) 

ml + k 2 = fi (35) 

—k 2 {m + l) = y (36) 

k 2 ml = 5 (37) 


From Eqs. (34) and (36), we obtain the magnitude of the imaginary roots that expresses the 
angular velocity to: 


to- = k 2 = - 
a 


(38) 


The angular velocity to as a function of S n and D u is given in the following equation: 


co = 


D dp {M, co, </>) • S pd (M) + S dp (M) • D pd - D dd {M, co, <j>) • S PP {M) - S dd {M) • D pp 1 1/2 


D dd {M,co, </>) +D 


pp 


Using Eqs. (37) and (38), Eq. (35) is written: 


y da. 


£ = - + -=► - ~P - +<5 = 0 

a y \a/ Va/ 

Substituting ^ as co 2 , the geometric constraint is obtained [9]: 
co 4 — ftco 2 + <5 = 0 


(39) 

(40) 


( 41 ) 
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Surface of Stable operation (SSO) 
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Fig. 2. Surface of stable operation (SSO) of RE-1000. 


The solution of the biquadratic equation using the coefficients S u and Du is 


co = 




1/2 


1/2 


(42) 


Furthermore, the parameters a, /?, y and 3 should be greater than zero, something that is usually 
satisfied. Eqs. (39) and (42) should have equal solutions, which is the acceptable angular velocity 
for stable operation of the FPSE. 


2.1.1. Case study I 

The previous analysis has been applied to a typical FPSE, the RE-1000 produced by Sunpower 
Inc. The diagram of angular velocity versus the displacer-piston amplitude phase angle and the 
working gas mass follows. 

As shown in Fig. 2, the effect of the phase angle on the angular velocity is much stronger than 
the effect of the working gas mass. This could be explained from the fact that the phase angle and 
angular velocity have effects on the same terms of the damping coefficients, while the working gas 
mass, except these, have effects on the stiffness coefficients. The diagram surface shows the stable 
operation area in a range from 0.5 to 3 gr working gas mass and from 50° to 75° phase angle. 


(i) In the area below the surface of stable operation (SSO), both the piston and displacer motions 
are unsteady and oscillate with increasing amplitudes. 

(ii) In the area above the SSO, the piston and displacer motions are oscillations whose amplitudes 
tend to zero. 


In Figs. 3 and 4, the stable operations of the engine are shown for a point on the SSO. 

In Figs. 5 and 6, the unstable operations of the engine are shown (case (b)—roots with only 
negative real parts—oscillations whose amplitudes tend to zero—stoppage). The greater the 
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Fig. 3. Stable operation. Piston velocity V p versus piston amplitude X p . 



Fig. 4. Stable operation. Displacer velocity Vd versus displacer amplitude X d . 


angular velocity is than the optimum (which is given by the SSO), the faster the displacer and 
piston amplitudes are reduced. 

In Figs. 7 and 8, the unstable operations of the engine are shown, (case (c)—roots with positive 
real parts are complex conjugates—there is an oscillation for each component, although they are 
divergent). Respectively, in this case, the less the angular velocity is than the optimum (which is 
given by the SSO), the faster the displacer and piston amplitudes increase. 


2.1.2. Case study II 

Typically, the operating values for the thermodynamic parameters of this engine, in which the 
working gas is helium, are 814 K heater temperature and 322 K cooler temperature. The above 
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Fig. 5. Unstable operation (stoppage). Piston velocity V p versus piston amplitude X p . 



Fig. 6. Unstable operation (stoppage). Displacer velocity Va versus displacer amplitude X A . 


analysis has been used while keeping these parameters constant, and the differential equations 
system has been solved parametrically for two modes. 

1. On the first mode, the independent parameter is the working gas mass, and the dependent 
parameters are the phase angle and the angular velocity. 

2. On the second mode, the independent parameter is the phase angle, and the dependent para¬ 
meters are the working gas mass and the angular velocity. 

The following Figs. 9-12 show how the working gas mass affects each operation parameter and 
other magnitudes, such as work done and indicated power output. 
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Fig. 7. Unstable operation (oscillation with increasing amplitude). Piston velocity V p versus piston amplitude X v . 
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Fig. 8. Unstable operation (oscillation with increasing amplitude). Displacer velocity V d versus displacer amplitude X d . 


In Fig. 9, the angular velocity versus the working gas mass is shown. On the solid line, there are 
the results of mode 1. The dot lines are the results of case study I and represent stable operation 
points for constant values of the phase angle. The points at which the solid line cuts the dot lines 
are the points at which mode 1 gives the same results as case study I. The solid line is represented 
approximately by the following equation: 

a = 158.2 -M 0 ' 55 (43) 

Fig. 10 shows how the piston and displacer amplitudes vary with the working gas mass. Also 
there is the ratio r of the amplitude X d to the amplitude X p . As shown, when the working gas mass 
is lower than 0.82 g, the piston and displacer amplitudes are constant, and the displacer amplitude 
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Fig. 9. Angular velocity versus working gas mass (for different phase angle values). 
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Fig. 10. Displacer amplitude, piston amplitude and ratio X A jX v versus working gas mass. 


is greater. For a working gas mass greater than 0.82 g, the displacer amplitude decreases and the 
piston amplitude increases. After the value of 1 g, it is clear that the piston amplitude becomes 
greater than that of the displacer. The ratio r is also constant in the beginning and then decreases 
according to the following polynomial equation: 

r = 0.81 + 0.909 • M - 1.01 • M 2 + 0.284M 3 (44) 

For calculations of the mean pressure of the engine, work done, heat transfer in cold and hot 
spaces and indicated power output, the isothermal analysis is used. Fig. 11 shows the mean 
pressure versus working gas mass and angular velocity. The mean pressure increases as the 
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Angular Velocity a> (rad/s) 

100 120 140 160 180 200 220 240 260 280 300 



0.5 1.0 1.5 2.0 2.5 3.0 

Working Gas Mass (g) 


Fig. 11. Mean pressure versus working gas mass and versus angular velocity. 



0.5 1.0 1.5 2.0 2.5 3.0 

Working Gas Mass (g) 


Fig. 12. Energy transfers and power versus working gas mass. 


working gas mass increases, while it decreases as the angular velocity decreases. The equations, 
which describe the mean pressure variation, follow: 

i?mean = 36.3 -M (45) 

Pmean = 1-13 + 0.038 • CO + 0.0012 • co 2 (46) 

The work done per cycle is equal to the sum of the heat transfer to the gas in the hot space per 
cycle and the heat transfer from the gas in the cold space per cycle. The work done and heat 
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transfers are shown in Fig. 12. On the right axis of the diagram, there is the indicated power 
output. The indicated power output increases as the working gas mass increases. The following 
equations represent the above magnitudes. 


W — 35.4 ■ M 

(47) 

Q e = 58.7 • M 

(48) 

Q c = —23.3 M 

(49) 

P = 0.546 • M + 0.367 • M 2 

(50) 


In contrast to mode 1, in mode 2, there are two acceptable values of each magnitude for each 
value of the phase angle. There is a minimum phase angle equal to the critical phase angle of 52.4° 
with critical angular velocity 140.3 rad/s, below which mode 2 does not find any solution for stable 
operation. For better comprehension, the results of this mode are separated into two groups. In 
Group 1 (Gl), the values of the magnitudes that correspond to the values of the angular velocity 
greater than the critical are included. In Group 2 (G2), the values of the magnitudes that cor¬ 
respond to the values of the angular velocity less than the critical are included. 

The following Figs. 13-16 show how the phase angle affects each operation parameter and 
other magnitudes, such as the work done and the indicated power output. 

In Fig. 13, the angular velocity and working gas mass versus the phase angle are shown. As 
shown in the diagram, the angular velocity in Gl increases at a greater rate than it decreases in 
G2. The equations for each group for angular velocity and working gas mass are the following. 

Gl : 

G2 : 



co = -4867.1 + 231.6 • </> - 3.6 • </> 2 + 0.02 • (51) 

to = 26808.1 - 1633.1 • </> + 37.4 • 4> 2 - 0.4 • <j) 3 + 1.4 x 10~ 3 • 4> 4 (52) 
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Fig. 14. Displacer amplitude, piston amplitude and ratio X A /X V versus angular velocity. 


Angular Velocity ca (rad/s) 

50 100 150 200 250 300 



Fig. 15. Mean pressure versus phase angle and versus angular velocity. 

Gl : M = -107.5 + 5.4 •</>- 9.1 x 10 2 • 0 2 — 5.2 x 1(T 4 • </> 3 (53) 


G2 : M = (—141.4 + 2.7 • ^>)~ 0 ' 243 ( 54) 

In Fig. 14, the piston and displacer amplitudes and ratio r of amplitude X d to amplitude X p are 
shown. Although the piston amplitude increases with angular velocity, the displacer amplitude 
has a maximum of about 1.45 cm at 103 rad/s. The ratio r is reduced continuously, beginning 
greater than 1.0 (displacer amplitude greater than piston amplitude), becoming equal to 1.0 at 
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Fig. 16. Energies (W, Q t , g c ) and power versus angular velocity. 


145 rad/s and continuing reducing. The two groups are separated with the dot line, as shown in 
the diagram. The relevant equation of the ratio r is the following: 

r = 2.12-0.011 -co + 3.02 x 1(T 5 • to 2 - 2.83 x 1(T 8 • to 3 (55) 

As shown in Fig. 15, the area of high values of angular velocity corresponds to high values of 
mean pressure. In Fig. 15, it is clear that the mean pressure increases as the angular velocity 
increases. The equations, which describe the mean pressure variation, are the following: 

Gl : Pmeaa = -3868.8 + 194.6 > - 3.3- (/> 2 + 0.02 • </> 3 (56) 


G2 : 


Pmean 


— 14061.1 + 279.3 • <fi 
1 - 14.8 • r/) + 0.29 • c/r 


(57) 


/’mean = 0-42 + 0.055 • CO + 0.001 • co 2 + 1.68 • to 3 (58) 

In Fig. 16, the work done per cycle, heat transfers in the hot and cold spaces per cycle and the 
indicated power output are shown. The following equations represent the above magnitudes: 

W = 2.65 + 0.008 • co + 0.001 • co 2 - 3.05 x 10~ 7 • co 3 (59) 

Q c = 4.34 + 0.013 • co + 0.002 • co 2 - 5.05 x 10~ 7 • co 3 (60) 

Q c = -1.74 - 0.005 • co - 8.6 x 10~ 4 • co 2 + 2 x 10~ 7 • co 3 (61) 

P = 0.084- 1.6 x lO^ 3 -co+1.68-to 2 +1.61-to 3 (62) 

In Table 1, the main parameters of the RE-1000 are shown, as it has been evaluated from Urieli 
and Berchowitz [9], Walker and Senft [3], the authors’ previous analysis [11], the two case studies 
of the present analysis and experiments [9]. All results seem to compare well with the phase angle, 
slightly better in case studies 1 and 2. The comparison with experimental results in the amplitude 
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Table 1 

Comparison of the operation parameters of RE-1000 


Parameters 

Urieli and 
Berch. 

Walker and 
Senft 

Previous 
analysis [11] 

Case study 1 

Case study 2 

Experiment 

Frequency (Hz) 

33.2 

32.9 

32.9 

30.0 

30.1 

30.0 

Phase angle (°) 

-57.9 

-55.1 

-55.1 

-54.8 

-55 

-42.5 

Amplitude ratio 

0.62 

0.63 

0.63 

0.87 

0.73 

1.06 


ratio is quite better in the present analysis, more specifically in case study 1. The same appears in 

the comparison of frequency, in which we can obtain the exact experimental value using the case 

study 1. 

3. Conclusions 

• The stiffness coefficients S i; are dependent only on the working gas mass M and not on the 
phase angle cj) or angular velocity to. For this reason, the constant order 8, of the characteristic 
equation of the differential equations system, is also dependent only on the working gas mass 
M. 

• The machine works at a cyclic steady state when the characteristic polynomial has two imagi¬ 
nary roots and two roots with negative real parts. The absolute values of these real parts must 
be very high to ensure a quick drop of the transient and, therefore, a quick arrival of the sta¬ 
tionary oscillations for the reciprocating parts. The farther an operation point is from the SSO, 
the faster the operation of the engine will be increased or decreased in each case. 

• The result of case study I is a surface of stable operation (SSO). In the area below the SSO, both 
the piston and displacer motions are unsteady and oscillate with increasing amplitudes. In the 
area above the SSO, the piston and displacer motions are oscillations whose amplitudes tend 
to zero. 

• As it seems in Fig. 8, the results of mode 1 of case study II are points of the SSO of case study I. 
The results of mode 2 are also points of the SSO. 

• In mode 1 of case study II, the piston and displacer amplitude are constant until 0.82 g working 
gas mass, with the piston amplitude greater. After this value, the piston amplitude is decreased, 
the displacer amplitude is increased, and as a result, the ratio r is decreased from 1.04 to 0.76. 
The ratio r becomes one when the working gas mass is approximately 0.9 g. 

• In contrast to mode 1, in mode 2, there are two acceptable values of each magnitude for each 
value of the phase angle. There is a minimum phase angle equal to the critical phase angle of 
52.4° with critical angular velocity 140.3 rad/s, below which mode 2 does not find any solution 
for stable operation. 

• The variation of mean pressure is greater in Group 1 of mode 2 than in Group 2, as is clearly 
seen in Fig. 14. The same behavior follows the work done per cycle and the power. 
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